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ABSTRACT 

Aims. This paper presents a new homogeneous catalogue of blue straggler stars (BSS) in Galactic open clusters. 

Methods. Photometric data for 216 clusters were collected from the literature and 2782 BSS candidates were extracted from 76 of them. 
Results. We found that the anticorrelation of BSS frequency vs. total magnitude identified in similar studies conducted on Galactic globular 
clusters extends to the open cluster regime: clusters with smaller total magnitude tend to have higher BSS frequencies. Moreover, a clear 
correlation between the BSS frequency (obtained normalising the total number of BSS either to the total cluster mass or, for the older clusters, 
to the total number of clump stars) and the age of the clusters was found. A simple model is developed here to try to explain this last and new 
result. The model allows us to ascertain the important effect played by mass loss in the evolution of open clusters. 
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1. Introduction 

Blue straggler stars (BSS) occupy a region above the turn-off 
(TO) point in the colour-magnitude diagram (CMD) of star 
clusters, where no stars are expected on the basis of standard 
stellar evolution, if we assume that all cluster stars are co- 
eval. A number of exotic explanations have been invoked to 
interpret the nature of these stars. Briefly, it has been sug- 
gested that these stars are the outcome of stellar collisions 
jBenz & Hillsl Il987t lLombardietal]ll996l) or of mass ex- 
change in close binary systems jMcCrealll964llEggen & Ibenl 
Il989t iMateo et alJ Il990l). Other possible s cenarios include 
second-generation stars feggen & Ibenl 1988I). accretion o f gas 
from the interstellar surrounding medium JWiihamsll96"4T) . and 
capture of field stars by a star cluster. 

Interestingly, BSS are ubiquitous, as they have been found 
in globular clusters (GCs), dwarf galaxies, open clusters (OCs), 
and in the Galactic field. Therefore, important constraints on 
their nature and possible formation mechanisms can be de- 
rived from a comparative study of the BSS properties in dif- 
ferent environments. In the past two decades, a huge amount of 
new, high-quality CMDs from CCD photometry from ground- 
based and s pace f acilities have become available. Recently, 
IPiotto et all J2004 .hereafter PDK04) published a new photo- 
metrically homogeneous catalogue of BSS extracted from the 



CMDs from HST observations of more than 1/3 of the known 
Galactic GCs JPiotto et"afll2002|) . The present paper aims to 
complement the PDK04 catalogue by extending the study of 
BSS properties to the different environment that is typical of 
OCs. 

The catalogue presented in this paper can be consid- 
ered as an important ext ension of the effort made by 
Ahumada & Lapasset ( 1995,hereafter AL95) to compile a cat- 
alogue of 969 BSS in 390 OCs, based on the photometry ob- 
tained with photographic plates and photoelectric detectors in 
the 60s, 70s, and 80s. The main conclusions drawn from their 
catalogue are that (1) BSS are present in clusters of all ages, 
(2) the fraction of BSS increases with the richness and age of a 
cluster, and (3) BSS show a remarkable degree of central con- 
centration. Moreover, the number of BSS over the number of 
main sequence stars seems constant up to an age of 400 Myr. 

The main reasons that led us to perform a new search of 
BSS in OCs are the following: (1) the AL95 catalogue was 
based on a highly inhomogeneous set of data, (2) a large num- 
ber of new, high-quality photometric studies of OC stars have 
appeared since AL95, and (3) the BSS selection criteria needed 
in our opinion to be revised, especially for young OCs. Most 
important, in order to estimate the relative number of BSS 
in different clusters, the total number of BSS must be nor- 
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malised to some reference stellar population in the same clus- 
ter. The choice made by AL95 to normalise the BSS counts to 
the main sequence stars creates some problems, especially in 
young clusters, where it is difficult to decide whether a star is a 
main sequence, a blue straggler, or a field star. 

The layout of this paper is as follows. In Sect. 2 we present 
the new database. In Sect. 3 we describe a few interesting prop- 
erties highlighted by the new catalogue. In Sect. 4 we develop a 
simple model in order to interpret the observed BSS frequency 
in OCs. Section 5 presents our conclusions. 

2. The new database 

We collected a database of B, V (in some cases I) band CCD 
photometry of 216 galactic OCs from the literature. The clus- 
ters in this new catalogue have ages ranging from 5.5 to 10 4 
Myr. From the entire sample we extracted a new catalogue of 
2782 BSS candidates. This new catalogue can be considered an 
updated version of the AL95 catalogue. Recent CCD photome- 
try was collected for 136 clusters out of the original 390 in the 
AL95 catalogue. Moreover, the database was extended with 80 
additional OCs for which CCD photometry is now available, 
but which were not included in the AL95 catalogue. 

Out of the whole sample of 216 open clusters, only 76 (the 
oldest ones) show BSS candidates. Fifty-nine of these clusters 
had sufficiently deep photometry to be considered suitable for 
the kind of analysis we were interested in, which involves esti- 
mating the total mass from the integrated magnitude (see Sect. 
13- li . Moreover, to avoid contamination effects from field stars, 
we decided to select for our analysis only those objects with a 
shorter distance from the cluster centre than the angular radius 
of the cluster given in the WEBDA archive 1 . After applying 
this selection, we were left with 774 BSS candidates. 

The list of these 59 clusters used in our analysis and their 
relevant parameters is found in Table ^ Column 2 gives the 
logarithm of the age, Col. 3 is the distance of the cluster in 
parsecs derived via isochrone fitting, Cols. 4 and 5 list the in- 
tegrated magnitudes (the first value is calculated based on the 
catalogue used for the present work, see Sect. 13.11 the second 
one from Lata et al., 2002), Col. 6 is the radius in arcmin from 
the WEBDA archive, Cols. 7 and 8 are the observed numbers 
of BSS and clump stars inside the cluster radius respectively. 
Column 9 lists the photometric bands used for this work (BV 
or VI). 

The entire catalogue (extracted from the old- 
est 76 clusters and without any selection in ra- 
dius) will be published in an electronic form at the 
CDS (http://cdsweb.u-strasbg.fri, and avail- 
able at the Padova Globular Cluster Group web pages 
( http : //dipastro . pd . astro . it/globular si. 

Table 2 (available on-line) lists all the 216 clusters analysed 
to extract the current catalogue. The following parameters are 
provided for each entry: right ascension and declination, red- 
dening, distance modulus, logarithm of age, apparent cluster 
radius, number of BSS candidates within the apparent radius 
and without limitation on the distance from the centre, number 

1 http : //ob swww . uni g e . ch/ webda/navi gat i on . html 



of clump stars, and reference to the adopted photometry (all the 
references are listed in the file refs.dat, also available on-line). 

Table 3 (available on-line) lists all the 2782 BSS candi- 
dates. Column 1 provides the name of the cluster where the 
candidate was found. The presence of an asterisk in column 2 
indicates the existence of a special note referring to this entry 
(all the notes are in the file notes.dat also available on-line). 
Cols. 3, 4, 5, and 6 give the identification number, the distance 
from the centre, the magnitude and the colour respectively for 
each candidate. 

The content of the electronic tables is described in details 
in the file ReadMe also available on-line. 

Although based on a smaller number of OCs than the one 
by AL95, the new catalogue contains a larger number of BSS 
and has the advantage of being much more homogeneous and 
therefore better suited to comparative analyses of BSS in dif- 
ferent environments. 

2.1. Selection criterion 

BSS occupy a region of the CMD usually bluer and brighter 
than the TO. To determ ine the TO position, we superimposed 
the lGirardi et al.l (12000) theoretical isochrones to the observed 
CMD. The isochrone fitting also provided us with an estimate 
of the apparent distance modulus and age for each cluster. 

To properly account for the photometric errors and the dif- 
ferential reddening, we defined a BSS selection region in the 
following way: 

- we defined a main sequence for binaries (V Bi , (B - V) Bi ) 
with equal mass components, shifting the theoretical zero- 
age main sequence (ZAMS) (V,-, (B - V),) 0.75 in magni- 
tude V toward brighter magnitudes 

V Bi = V t - 0.75 
(B - V) Bi = (B- V)i 

- we then defined the two borders (I for left and r for right) 
of the selection region 

Vu = Vi + AV 
(B - V)u = (B- n - MB - V) 

V ri = V Bi - AV 
(B - V) ri - (B - V) Bi + A(B - V) 

where AV and A(B - V) were set to 0. 15 magnitudes to take 
photometric errors and differential reddening into account. 

Figure^shows an example of the selection region on the CMD 
of NGC 7789. The same procedure has been applied to the 
CMDs in (V, V-I). 

3. Analysis 

This section presents our analysis of the selected BSS cata- 
logue. The estimate of the main parameters used in this study 
is explained in details and the most relevant correlations among 
the observed quantities are described. 
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Fig. 1. Colour-magnitude diagram of the open cluster 
NGC 7789. Solid lines: theoretical ZAMS for single stars and 
for binaries with equal-mass components. Dashed lines: same 
sequences shifted by 0. 15 magnitudes in colour and magnitude 
to take photometric errors and differential reddening into ac- 
count. Filled circles: BSS candidates. 



3.1. Integrated magnitudes 

For each cluster we computed the integrated magnitude and the 
total mass, considering as cluster members only the stars falling 
inside a certain angular radius R from the centre (the adopted 
values of R are those listed on WEBDA and in Tabled 2 . 

We calculated the total magnitude M v 5 as the sum of the 
luminosity of all the stars with absolute magnitude My < 5 and 
angular distance from the centre smaller than R. The adopted 
limiting absolute magnitude is a compromise between the ne- 
cessity of reaching faint magnitudes, in order to include as 
many cluster stars as possible, and of having enough clusters 
with deep enough available photometry to guarantee a homo- 
geneous determination of the total cluster magnitude. 

Recently. lLata et aD ll2002l hereafter LP02) published a cat- 
alogue of integrated magnitudes for a large sample of open 
clusters. Thirty-eight out of 59 clusters in our sample were 
listed in their catalogue. Figure |2 shows the comparison of 
LP02 values with respect to our values for the 38 clusters in 
common. There is overall agreement, although our integrated 
magnitudes are systematically slightly brighter, probably due 
to contamination. 



Fig. 2. Comparison between our integrated luminosity esti- 
mates and those published by LP02. The solid line merely 



shows the relation M v 5 



M 



LP02 



2 Note that the lack of unresolved stars in our total magnitude es- 
timate is partly compensated for by the likely presence of some field 
stars left in the bright star catalogue. 



3.2. Observed number of BSS 

The upper and lower panels of Fig.[3]show the trend in the ob- 
served number of BSS stars (Ness) in each cluster vs M v 5 and 
My P02 , respectively. This figure is an extension t o lower lumi- 
nosity and into the OC mass regime of Fig. 1 in iDavies et alJ 
The number of BSS is generally less in OCs than in 
GCs. Though with larger dispersion, the number of OC BSS 
versus total magnitude seems to lie on the prolongation of the 
trend indicated by the l ower-luminos ity GCs and by the mod- 
els developed bv lDavies et alJ (J2004) (Fig.EJ. We would like 
to stress here that the model developed bv lDavies et alJ (|2004) 
is not applicable to the OC environment. It is shown here only 
as a useful comparison. The larger dispersion of the OC data 
with respect to the GC ones is due to the spread in age. Older 
clusters have a larger number of BSS. This effect could not be 
observed in the GC case, as the spread in age is much less sig- 
nificative. The effects of age on the BSS population in OCs will 
be discussed further in the next sections. 



3.3. Specific BSS frequencies 

In general, we can expect that more massive clusters contain, 
on average, more stars in any evolutionary branch, including 
BSS. Therefore, to properly study the BSS properties in differ- 
ent clusters, it is common practice to normalise the number of 
BSS to some global cluster parameter or to the total number of 
stars in some well-defined region of the CMD. 

We considered two different normalisations in order to es- 
timate the relative frequency of BSS in each cluster: 
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Fig. 3. Correlation between log Nbss and the integrated mag- 
nitude. Circles represent OC data, open symbols are used to 
mark clusters were only 1 or 2 BSS candidate were found. The 
GC data from PDK04 are plotted using triangles. The solid line 
shows the theoretical prediction by Davies et al., 2004. 



In the following, we discuss the relevant correlations that 
were found between the derived BSS frequency and other clus- 
ter parameters. 

- Normalisation with total mass 

Figure|4]shows the logarithm of the specific frequency de- 
rived by normalising the total number of BSS to the cluster 
total mass (Nbss /Mm) plotted as a function of the loga- 
rithm of age. There is a clear correlation between these two 
quantities: older clusters have a larger relative number of 
BSS. A marginal correlation between these two quantities 
has already been found by AL95; however, we believe that 
this is now much more evident with the new catalogue. 
PDK04 show that the relative number of BSS in GCs an- 
ticorrelates with the cluster total luminosity, in the sense 
that less luminous (less massive) GCs have a larger rela- 
tive number of BSS. Recently, this anticorrelatio n has been 
further extended to very low-luminosity GCs by Sandauist 
Figure [5j shows the specific frequency of BSS vs 
the integrated magnitude. The anticorrelation with the to- 
tal magnitude continues down to M v 5 < -5: less massive 
clusters have, on average, a higher BSS frequency. In the 
smallest clusters of our sample, the BSS relative frequency 
is about two orders of magnitude higher than in the most 
massive GCs ofPDK04. 




8.5 
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Fig. 4. The figure shows the trend in log(A^BSs/A1tot) vs the log- 
arithm of age. Open symbols are used to mark clusters with 
only 1 or 2 BSS candidates. 



- A^BSs/Altot : number of BSS normalised to the total mass 
(expressed in solar masses Ale), derived from the total lu- 
minosity, assuming Attot/Xtot = 1-5 M /£o 

- Nbss/N c i : number of BSS normalised to the number of 
clump stars (i.e., stars on the horizontal branch, when this 
feature was clearly identified in the CMD). 



- Normalisation with clump stars 

The selection of clump stars in our CMDs was often very 
uncertain. The clump region in OCs is much less clearly 
defined and much less populated than the horizontal branch 
region in GCs. Only 31 out of the 59 clusters show a clear 
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Fig. 5. The figure shows the trend in log(A^BSs/A1tot) vs the in- 
tegrated magnitude. Open symbols are used to mark clusters 
with only 1 or 2 BSS candidates. 



clump population. These are the oldest clusters in our sam- 
ple (all older than ~ 10 9 years). The number of clump stars 
in every single cluster of our database is always smaller 
than 10. Nevertheless we do not expect this number to be 
contaminated by field stars, given the very small area occu- 
pied by clump stars in OCs. 

We then divided the subsample of 31 clusters into two 
groups, one containing 17 clusters where the selection of 
clump stars was straightforward and reliable, and another 
containing 14 clusters where the selection was uncertain 
mainly due to the small number of clump stars. 
Figure[6]shows the trend of log(A^BSs/^V c i) vs M y 5 , and the 
GC data from PDK04, who calculated the BSS frequency 
by normalising the total number of BSS to the total number 
of horizontal branch stars, which is analogous to our nor- 
malisation to the clump stars. Noteworthy, and even more 
clearly so than in Fig. [3] the trend in relative number of 
BSS vs integrated magnitude already se en for the glo bular 
cluster s also extends to the OCs. Finally. IPreston & Snedenl 
(2000) have estimated a ratio Nbss /Nub ~ 4 among field 
stars that agrees with the extrapolation of the correlation in 
Fig-EJat fainter integrated magnitudes. 

4. Modelling the expected number of BSS in OCs 

iDavies et alJ {2004) interpret the anticorrelation between the 
frequency of BSS and the total cluster mass in terms of the 
evolution of primordial binaries, which is affected by the stel- 
lar encounters. However, the previous section also shows a 
clear correlation between the number of BSS normalised to 
the total mass of the cluster and the age of the cluster itself, 




Mf 



Fig. 6. The figure shows the trend in log(A^BSS /N c \) vs the inte- 
grated magnitude. Circles represent OC data and open symbols 
have been used to mark clusters where the selection of clump 
stars was particularly uncertain. The GC data from PDK04 are 
plotted using triangles. 



with younger clusters having a lower relative number of BSS 
with respect to the older ones. The GCs in the PDK04 cat- 
alog ue were all coeva l, with ages ranging from 8 to 11 Gyr 
(iDe Aneeli et 311200 5*). The OCs have a much wider spread in 
age, from 10 8 to 10 10 yr. This property offers some insight into 
the possible connection between cluster evolution and BSS for- 
mation. 

4.1. Model description 

In order to further investigate this result and in an attempt to 
interpret it, we have developed a model following some simple 
prescriptions. We first assume that in OCs the number of BSS 
formed through direct collision (dynamical BSS) is negligible 
compared to the number of those formed through the evolu- 
tion of primordial binaries (primordial BSS). This assumption 
is reasonable when considering the following facts: 

- In the first instance, one could consider the estimated num- 
ber of encounters involving single stars. The time-scale for 
a single star to under go an encounter with anoth er star may 
be approximated by (Binnev & Tremaine 1987) 

where n is the typical number density (calculated as the to- 
tal number of stars divided by the volume in pc 3 ), and M 
the typical combined mass of the two colliding stars (as- 
sumed to be of the order of 2 -f 3 M Q ). The value of V ra can 
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be calculated for each cluster as V2GAt c iuster /^cluster- The 
ratio between the age of each cluster (t c £) and T enc 



(2) 



should give a reasonable estimate of the number of encoun- 
ters that could produce BSS. The ratio N snc /Miot is negligi- 
ble compared to the number of observed BSS (the expected 
number being between 1 and 2 orders of magnitude lower 
than the observed number of BSS). Moreover, this is likely 
to overestimate the number of collisional BSS because 1) 
not all the dynamical interactions will produce BSS, 2) only 
BSS formed in the past m 1 Gyr (approximate estimate of 
the BSS lifetime) will be visible today. On the other hand, 
this estimate does not take dynamical interactions involv- 
ing binaries into account. These are expected to be more 
frequent due to the larger cross section. However, even in 
this case, we do not expect all the encounters to produce 
BSS. 

- The fraction of binaries undergoing a dynamical interaction 
can be estimated using the former formula but adopting a 
value for R m j n typical of the size of a binary. This will ob- 
viously depend on the properties of the primordial binary 
population, and will be higher for low-velocity dispersion 
clusters (as the velocity dispersion decreases, wider bina- 
ries qualify as "hard"). Obs ervational evidence (for a dis- 
cussion see ISandauistll2005l) shows that at least » 25% of 
the BSS in M 67 are dynamically formed. However, assum- 
ing a value of the order of a 100 AU for R m ; n and using a 
realistic estimate of the cluster velocity dispersion in Eq.^ 
one finds that the mean expected fraction of binaries experi- 
encing dynamical interactions for the clusters in our sample 
is more than 10 times lower than the one in M 67. Thus we 
expect the fraction of dynamical BSS in most of our clus- 
ters to be of the order of a few percent, assuming an equal 
binary fraction for all clusters. More careful considerations 
should also take the age of the clusters into account. Our 
sample spans an age range between 10 8 and 10 ll) yr, and it 
is not at all clear how age influences the collision rate. The 
widest binaries (those with higher cross-section) are likely 
to interact earlier on, and their BSS products might not be 
visible today. 

Our simple model builds up a cluster with both a sin- 
gle star and a binary population. All the components are 
drawn randomly from the same IMF (we adopted the IMF of 
Egg leton et alJl989T) . The masses are generated using the equa- 
tion 

19jc 

M = (1-^^0.032(1 -^ ^' (3) 

with a series of random numbers x. The binary fraction is as- 
sumed to be of the order of 0.5. 

Binary systems will form BSS only if one of the two fol- 
lowing criteria is verified: 

1 . either the mass ratio between the secondary and the primary 
(q - Mz/Mi) is higher than a certain q MT (and obviously 
lower than 1); in this case the mass transfer between the 
primary and the secondary is expected to be stable and to 



lead to the formation of a blue straggler star with an evolved 
companion; 

2. or the mass ratio between the secondary and the primary 
is lower than ^mt but higher than a certain qcE, in which 
case the mass transfer will probably be unstable and lead to 
a merger of the two components of the binary system, this 
time leading to the formation of a single blue straggler star. 

If q < qcE, we believe that a common envelope will form and 
no BSS will originate from such systems. 

The choice of the values for <?mt and qcE was taken in con- 
sideration of the following facts. Conservative mass transfer 
from a more-mass ive donor to a l ess-massive receiving star 
is often unstable jFrank et alJl2002l) . However, if the mass ra- 
tio is close to unity, mass transfer will reverse the mass ratio 
(i.e. the donor will quickly become the less-massive star) and 
mass transfer may proceed in a stable fashion. There is evi- 
dence that this must be the case for at least some BSS, since 
some field BSS ha ve been found in w ide binaries presumably 
with white dwarfs JPreston & Sneder]l2000l) . As a reasonable 
example here, we take this critical mass ratio to be ^mt = 0.85. 

In systems with a mass ratio below this value, mass trans- 
fer will be unstable, resulting in the merger of the two stars and 
the formation of a blue straggler, if the mass exceeds the current 
turnoff. In systems with a very low mass-ratio, the less-massive 
star will be much denser than the primary. Rather than merge, 
the two stars will form a common envelope system where gas 
from the more-massive (and less dense) envelope will smother 
the less-massive (and denser) star, forming a common enve- 
lope around it and the core of the more-massive star. In this 
second case, the common envelope will be ejected as the core 
and the low-mass star spiral together, and no blue straggler 
will be formed. Here we assume this occurs for q < 0.4 (i.e. 
qcE = 0.4). 

For each binary, leading to a blue straggler, three different 
times are defined: 

1. toN : the time at which the primary evolves off the main 
sequence, and thus presumably starts the mass transfer on 
the secondary or the merger between the two companions. 
The main sequence lifetimes were calculated by applying 
the analytic formulae for stell ar evolution as a fu nction of 
mass and metallicity given bv lHurlev et alJ J2000t) . assum- 
ing a typical metallicity close to the solar one. 

2. tess : the lifetime of the formed star. This is calculated in 
a different way depending on the formation process: (1) for 
mass transfers, this is the main sequence lifetime (t) of a 
star of mass Mi + 5/6 Mi as the expected transferred mass 
is about 5/6 of the mass of the primary; (2) for mergers, 
this will have to be calculated considering the fraction of 
the core mass left available in the two merging stars. In 
particular, if Mi and M2 are the starting total masses of 
the two components, then at the time T\ when the primary 
evolves off the main sequence, we will have that the core 
mass of the primary 

lb 



M c 



= aM\ 



(where a = 0.125 and b - 1 .4. Ivan den Heuvellll994 has 
already been transformed into helium (Alcore.i = Mhci)- 
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At the same time, assuming that the hydrogen burning pro- 
ceeds linearly with time during the evolution on the main 
sequence, the secondary will have burned 

AW = aM b 2 TilT 2 = Mcore,2Tl/T 2 . 

The core mass of the merger product will be At core ,i+2 = 
a(Mi + Mif, but not all of it will be hydrogen. In fact 
the hydrogen core mass will be Mh,i+2 = Af C0le ,i+2 - 
AtHe.i ~ AlHe,2- The main sequence lifetime will then be 

T l+2A1 H ,l+2/A1 core ,l+2- 

3- toFF = toN + tBSS : the expected time at which the formed 
star evolves off the main sequence, leaving the CMD area 
where BSS are selected. 

Finally, the model loops over all the time range of interest 
(10 8 < t < 10 10 Gyr) and for each time (i.e. age) counts all the 
BSS having toN < t < toFF- At each time t, the total mass of 
the cluster is estimated as the sum of all the unevolved single 
stars and binaries (we consider only those binaries in which the 
primary has not yet evolved off the main sequence, t\ < t). The 
total mass will then be used to normalise the total number of 
BSS to allow a comparison with the observed frequencies. 

4.2. Interpreting the empirical results 

The resulting expected time evolution of the relative number 
of BSS normalised to the total mass is only slightly depen- 
dent on the assumed values for ^mt and qcE- Changing the 
values of these parameters shifts the expected trend toward 
higher or lower frequencies but does not influence the evo- 
lution rate significantly. Clearly many uncertainties must be 
taken into account when proceeding with this kind of anal- 
ysis. Observationally, the uncertainty on the mass-luminosity 
ratio and thus on the total mass clearly affects the level of the 
trend. Theoretically, the adopted values for the parameters qmt 
and qcE can also shift the expected trend towards higher or 
lower frequencies, and our values, although reasonable, have 
not been precisely determined in previous works. For these rea- 
sons, we concentrate in the following more on the general trend 
of the correlation than on the level of the relative numbers of 
BSS. Figure [71 shows the expected trend assuming qmj = 0.85 
and qcE = 0.4, with the corresponding observed quantities. As 
shown in the plot, the expected trend of BSS frequency with 
age is flatter than the observed one. Older OCs show a higher 
number of BSS per unit mass than expected. 

However, one important aspect has not been taken into con- 
sideration by our model. Mass loss plays a very important role 
in OCs, so we should account for it when defining the total nor- 
malising mass. Since they are more massive than normal stars, 
BSS sink toward the centre quickly after formation, which im- 
plies that the BSS population should be less affected by mass 
loss than the normal stellar population. At increasing ages, OCs 
tend to lose more and more mass. Therefore the normalising 
factor becomes smaller and smaller, steepening the relation be- 
tween normalised number of BSS and cluster age in Fig. Q 
Mass loss also depends on the cluster's total mass, which can 
partly account for the wide spread of the observed quantities, 
together with the large uncertainties in the mass estimates. 




Fig. 7. Comparison between expected log(A^BSs/A4tot) and ob- 
served frequencies of BSS in OC. The solid line represents the 
result of our model. Observed data are shown with filled and 
open circles, open circles marking clusters with only 1 or 2 
BSS candidates. 



Recently, a couple of papers h ave addressed the is- 
sue of mass loss in stell ar clusters. [ Tanikawa & Fukushied 
J2005t hereafter TF05) and lLamers et alJ J2005t hereafter L051 
present Nbody simulation results, predicting mass loss 
timescales of star clusters. We applied the results of both works 
to investigating the effect of mass loss on our predictions. 

Figure 2 in TF05 reproduces the mass loss as a function of 
time. The curves from the first panel at the top and on the left 
are the most appropriate for our analysis, having lower initial 
concentration and higher resolution. Time was converted from 
Nbody units to Myr by adopting the typical cluster properties 
in 5 different age intervals. This gave us the typical mass loss 
in each interval, which we then applied to the results of our 
model, assuming that BSS are not affected by mass loss. We 
repeated this procedure for all the curves in the first panel of 
Fig. 2 of TF05 (2k, 8k, 32k, 131k, where k stands for 1000 
particles), and Fig.[8]shows the different results. The estimates 
obtained using the 2k and 8k Nbody models are probably the 
more realistic ones, all our clusters having less than 2100 stars. 

The observed points nicely occupy the region between the 
lines, representing the two most extreme scenarios (no mass 
loss, and large mass loss affecting only normal stars), with only 
a few points falling outside this region. This may suggest that 
BSS are indeed affected by mass loss as well, although not as 
much as normal stars. This preliminary analysis confirms that 
our simple model with approximated mass-loss estimates can 
reproduce the distribution of the observed points quite well. 

To confirm these conclusions, we tried to estimate the mass 
loss with a different approach. L05 provide a set of analytical 
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Fig. 8. Comparison between expected log^BSs/Altot) and ob- 
served frequencies of BSS in OC. The solid line represents the 
result of our model. Observed data are shown with filled and 
open circles, open circles marking clusters with only 1 or 2 
BSS candidates. Dotted, dashed, dot-dashed and long-dashed 
curves reproduce respectively the expected values obtained tak- 
ing into account mass loss following the results of Tanikawa & 
Fukushige (2005) with 2k, 8k, 32k, and 131kNbody models. 



Fig. 9. Comparison between expected log(A^BSs/Wl to t) and ob- 
served frequencies of BSS in OC. The solid line represents the 
result of our model. Observed data are shown with filled and 
open circles, open circles marking clusters with only 1 or 2 
BSS candidates. Asterisks are the observed BSS frequencies af- 
ter correcting for mass loss due to tidal effects using the recipe 
suggested by Lamers et al. (2005) as explained in the text. 



formulae to estimate the mass loss due to both stellar evolu- 
tion and tidal effects. Stellar evolution is taken into account by 
our model. Due to the analytical form of the recipe provided 
by L05, it is convenient to eliminate the tidal effects from the 
observed data instead of adding them to the expected quanti- 
ties. For this reason, we use Eq. 7 of the L05 paper to estimate 
the initial mass of the clusters in our catalogue starting from 
the present observed mass (note that this equation takes into 
account the mass loss due to both stellar evolution and tidal 
effects). Then using Eqs. 3 and 2 of L05, we correct the calcu- 
lated initial mass taking the mass loss due only to stellar evolu- 
tion into account. The obtained mass estimates will be affected 
by mass loss due to stellar evolution but not by mass loss due 
to tidal effects, i.e. will be comparable with the masses calcu- 
lated in our model (which does not account for tidal effects but 
includes stellar evolution). 

Figure[9]shows the observed values, the corresponding val- 
ues corrected for mass loss due to tidal effects applied to normal 
stars only and our model predictions. The corrected values rep- 
resent the case where the mass loss is the most evident. If BSS 
were also affected by mass loss, then the observed frequency 
would lie between the observed values and the corrected ones. 
The trend predicted by the model falls in between the observed 
quantities and those corrected for the maximum tidal effect. 
Clearly the results obtained with both approaches are consis- 
tent with each other thus suggesting that tidal effects play an 



important role in OC and that they probably affect the number 
of observed BSS, as well as the total mass. 

In this respect we wish to make one last consideration. The 
fraction of binaries undergoing a dynamical interaction can be 
estimated using Eq. ^ as explained in Sect. 4. 1 . This fraction 
can be estimated for each cluster in our sample. For older clus- 
ters the expected values are higher than those obtained for the 
younger ones. All but one of our clusters have an estimate 
for the fraction of binaries undergoing a dynamical interaction 
lower than the one calculated for M 67, the mean value being 
about one order of magnitude lower. Remembering that at least 
25% of the BSS population of M 67 appears to be dynamically 
formed, this could imply that several of our clusters could show 
a non negligible dynamical population of BSS. This would also 
cause the frequency of BSS in older clusters to be higher than 
the predictions of our model. 

Finally, the normalisation with the total mass adopted for 
the observed number of BSS, is sensitive to the choice of the 
mass-luminosity ratio. The dependence of A1 to t/-Ctot vs a g e is 
affected by both evolutionary and dynamical effects. Both stel- 
lar evolution and mass loss due to tidal interactions are ex- 
pected to cause the mass-luminosity ratio in old clusters to be 
higher than in young clusters. This would thus bring the ob- 
served trend in better agreement with the results of our model. 

As already noted in the previous section, there is a clear 
continuity betw een open and globular clusters. This is clearly 
shown in Fig. [3] D avies et alJ J2004I) use a model similar to the 
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one implemented here to predict the trend of the total number 
of BSS with respect to the total magnitude. 

Although in this paper we concentrate on the evolution of 
the BSS population with cluster age, we have also checked 
that the present m odel is consistent with the predictions of 
iDavies et alJ d2004h at the low mass limit of our OCs. At ages 
of the order of the Galactic GC ages (« 10 Gyr), the present 
model predicts logA^BSS - 3.2. From the total mass, predicted 
by the model for clusters with this age, we can estimate the to- 
tal magnitude assuming a mass-to-light ratio of the order of 3, 
as in lDavies et all J2004I) . We obtain M v « -9. At this magni- 
tude, their model predicts log jVbss ~ 2.1, including the effects 
of both primordial binary evoluti on and stella r collisions. The 
dashed curve in Fig. 6 of IDavies et all {2004) shows the con- 
tribution to the total number of BSS from primordial system. 
This curve would be a straight line if there were no collisions; 
in fact, the decline for magnitudes bri ghter than ~7-4 is d ue to 
the increasing 2V e nc. As discussed in Davi es et alJ d2004l) . the 
effect of the encounters is to fasten the evolution of primordial 
binaries: the more encounters we have the smaller the number 
of binaries able to evolve in a BSS after 10 Gyr or so. But this 
does not apply to the case of open clusters. As shown above, we 
do not expect encounters to play an important role in this envi- 
ronment, so that in this case we would expect the dashed curve 
to be linearly increasing with My, reaching logA^ess ~ 3.2 at 

My * -9. 

One further constraint must be satisfied by our model, 
i.e. the observed freq uency of BSS among field stars. 
IPreston & Snedenl ( 120001) estimate s BS s = N B ss/Nhb = 4 in 
the field. In order to compare our results with this value, we 
estimated the number of horizontal branch (HB) stars in our 
simulated cluster s. We again adopted the analytic formulae for 
stellar evolution faurlev et alJ koOO) to evolve our simulated 
stars and to count the number of stars during the HB phase 
at each time step. Assuming then the mu ch lower metallic- 
ity (Z = 0.0004) typical of field stars as in lPreston & Snedenl 
(2000), we recomputed our expected frequency of BSS with re- 
spect to the number of HB stars . For ages of the order of 1 2 -f 1 3 
Gyr, we obtained sbss ~ 2.5 4- 3, slightly lo wer but of the same 
order o f magnitude as the results obtained bv lPreston & Sned"en1 
(2000). The difference can be due to the different model, stellar 
evolution prescriptions, and input quantities. Obviously, when 
estimating the frequency of BSS in the field, mass loss should 
no longer play an important role, and thus the correct estimate 
should be given by our model without any correction for mass 
loss. 
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Table 1. List of open clusters and relevant parameters. 



ID 


logf 


Dist. (pc) 


Mf 




R 


A^BSS 




Band 


ArpMadore 2 


9.40 


13320 


-3.48 




1.00 


2 


3 


BV 


Berkeley 10 


8.80 


2294 


-5.21 




5.00 


7 




BV 


Berkeley 12 


9.50 


3289 


-2.67 




2.00 


14 




BV 


Berkeley 14 


9.20 


5495 


-6.54 


-4.07 


2.50 


40 


5 


BV 


Berkeley 20 


9.78 


8387 


-2.06 




1.00 


4 




BV 


Berkeley 21 


9.34 


5012 


-4.91 


-3.41 


2.50 


17 




BV 


Berkeley 22 


9.35 


6677 


-3.00 


-3.27 


0.50 


3 


2 


BV 


Berkeley 23 


8.90 


6918 


-4.91 




2.00 


3 




BV 


Berkeley 29 


9.50 


15180 


-3.77 


-4.64 


3.00 


3 


2 


BV 


Berkeley 30 


8.48 


4786 


-4.50 


-5.11 


1.50 


1 




BV 


Berkeley 3 1 


9.31 


7185 


-3.81 


-2.57 


2.50 


13 


4 


BV 


Berkeley 32 


9.53 


3496 


-4.00 


-2.97 


3.00 


19 


4 


BV 


Berkeley 33 


8.70 


3855 


-4.47 


-4.61 


2.00 


1 




BV 


Berkeley 39 


9.78 


3802 


-3.69 


-4.28 


3.50 


18 


5 


BV 


Berkeley 81 


9.00 


3006 


-4.74 


-4.70 


2.50 


9 




BV 


Berkeley 99 


9.50 


4943 


-4.36 


-3.66 


5.00 


6 




BV 


Collinder 74 


9.11 


2514 


-3.85 


-2.88 


2.50 


5 




BV 


Collinder 110 


9.11 


2415 


-5.07 




9.00 


18 


7 


BV 


Collinder 261 


9.85 


2518 


-3.92 


-2.87 


4.50 


35 


6 


BV 


IC 1311 


9.20 


6289 


-4.02 


-3.90 


2.50 


3 




BV 


IC 4651 


9.05 


882 


-1.20 


-2.89 


5.00 


1 




BV 


King 2 


9.78 


6003 


-3.41 




2.50 


20 


5 


BV 


King 1 1 


9.55 


2500 


-3.59 




2.50 


10 




BV 


Melotte 66 


9.60 


4018 


-6.40 




7.00 


29 


9 


VI 


Melotte 105 


8.32 


2210 


-4.85 


-4.27 


2.50 


1 




BV 


NGC 188 


9.63 


2047 


-3.91 


-2.86 


8.50 


30 




BV 


NGC 1193 


9.90 


5207 


-3.02 


-3.00 


1.50 


11 




BV 


NGC 1245 


8.90 


3012 


-4.50 


-3.60 


4.50 


7 


5 


BV 


NGC 1907 


8.57 


1558 


-4.26 


-6.19 


3.50 


2 




BV 


NGC 2112 


9.30 


850 


-3.04 


-1.88 


9.00 


6 




BV 


NGC 2141 


9.40 


3802 


-4.01 




5.00 


8 


6 


BV 


NGC 2158 


9.30 


3899 


-5.55 


-2.81 


2.50 


28 


8 


BV 


NGC 2194 


8.75 


2805 


-5.86 




4.50 


6 




BV 


NGC 2204 


9.20 


3981 


-5.47 


-4.65 


5.00 


6 


6 


BV 


NGC 2243 


9.45 


3976 


-3.26 


-2.67 


2.50 


7 


3 


BV 


NGC 2266 


9.00 


3148 


-3.99 


-3.83 


2.50 


2 




BV 


NGC 2420 


9.30 


2449 


-2.90 


-3.44 


2.50 


3 


1 


BV 


NGC 2477 


9.00 


1259 


-5.11 


-5.70 


7.50 


6 




BV 


NGC 2506 


9.25 


3089 


-4.51 


-4.31 


6.00 


8 


5 


BV 


NGC 2627 


9.15 


1858 


-3.59 


-4.84 


4.00 


7 


3 


VI 


NGC 2660 


9.03 


2328 


-3.79 


-4.08 


1.75 


5 


6 


BV 


NGC 2682 


9.50 


865 


-2.93 


-3.16 


12.50 


14 


3 


BV 


NGC 5999 


8.60 


2046 


-3.53 


-3.85 


1.50 


2 




BV 


NGC 6005 


8.95 


2018 


-4.50 


-3.66 


2.50 


6 


6 


BV 


NGC 6067 


8.10 


2241 


-6.66 


-6.18 


7.00 


1 




BV 


NGC 6253 


9.48 


1727 


-3.09 


-2.68 


2.00 


10 


4 


BV 


NGC 6791 


9.95 


3695 


-5.10 


-4.14 


5.00 


119 




BV 


NGC 6819 


9.38 


2349 


-4.10 




2.50 


5 


4 


BV 


NGC 6939 


9.06 


1862 


-4.13 




5.00 


9 


7 


BV 


NGC 7044 


9.10 


2989 


-4.59 


-3.93 


3.00 


9 




VI 


NGC 7789 


9.23 


1888 


-5.19 




12.50 


17 




BV 


Pismis 2 


9.06 


3311 


-5.13 




2.00 


17 


7 


BV 


Pismis 3 


9.03 


1361 


-4.11 




2.50 


12 


6 


BV 


Pismis 18 


9.00 


1941 


-3.89 




2.00 


2 


4 


BV 



Table 1. continued. 
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ID 


log* 


Dist. (pc) 


M< 5 




R 


A^BSS 




Band 


Saurer 2 


9.30 


6274 


-5.48 




2.00 


54 


6 


VI 


Tombaugh 1 


9.20 


3006 


-4.07 




2.50 


3 




VI 


Tombaugh 2 


9.20 


9419 


-5.14 




1.50 


7 


6 


BV 


Tombaugh 5 


8.30 


1754 


-6.06 




7.00 


3 




BV 


Trumpler 5 


9.50 


2995 


-5.88 


-5.06 


7.00 


60 


14 


BV 



